iV-dimensional Vaidya metric with cosmological constant in double-null coordinates 
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A recently proposed approach to the construction of the Vaidya metric in double-null coordinates 
for generic mass functions is extended to the n-dimensional (n > 2) case and to allow the inclusion 
of a cosmological constant. The approach is based on a qualitative study of the null-geodesics, al- 
lowing the description of light-cones and revealing many features of the underlying causal structure. 
Possible applications are illustrated by explicit examples. Some new exact solutions are also pre- 
sented and discussed. The results presented here can simplify considerably the study of spherically 
symmetric gravitational collapse and mass accretion in arbitrary dimensions. 
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I. INTRODUCTION 

The Vaidya metric [l| is a solution of Einstein's equa- 
tions for a spherically symmetric body with a unidirec- 
tional radial null-fluid. It has been used in the analysis 
of spherically symmetric collapse and the formation of 
naked singularities for many years (For references, see 
the extensive list of Q and also @). It is also known 
that Vaidya metric can be obtained from the Tolman 
metric by taking appropriate limits in the self-similar 
case [i]]- This result has shed some light on the nature of 
the so-called shell- focusing singularities [j|, as discussed 
in details in 0, H, 0, @- The Vaidya metric has also 
proved to be useful in the study of Hawking radiation 
and the process of black-hole evaporation [c| [111, H3 |. 
in the stochastic gravity program [13|, and, more recently, 
in the quasinormal modes analysis of varying mass black 
holes (3I3- 

The n-dimensional Vaidya metric is required in many 
physically relevant situations. The study of the gravita- 
tional collapse in higher dimensional spacetimes[16||. for 
instance, has contributed to the elucidation of the forma- 
tion, nature and eventual visibility of singularities. This 
last topic belongs to the realm of the celebrated Pen- 
rose's Cosmic Censorship Conjecture, see, e.g., [|[ for 
references. Higher dimensional varying mass black holes 
arc also central protagonists in the new phenomenological 
models with extra dimensions 17| . These objects might 
be obtained in the LHC at CERN[lc| and, once pro- 
duced, they are expected to decay driven by the emission 
of Hawing radiation. The analysis of their quasinormal 
modes can help the understanding of the dynamics of the 
evaporation process, and could even lead to some obser- 
vational signs [HI]. 

The inclusion of a cosmological constant A in the n- 
dimensional Vaidya metric is mainly motivated by the 
recent intensive activities in the AdS/CFT and dS/CFT 



conjectures, see [19j for some references. A cosmological 
constant is also necessary to allow the discussion of the 
n = 3 case in the same framework of the higher dimen- 
sional cases. 

The n-dimensional Vaidya metric was first discussed in 
[2(il ]. It can be easily cast in n-dimensional (n > 3, the 
three dimensional case will be discussed latter) radiation 
coordinates (w,r,8i, . . . , #„_2)[lll: 



ds = — I 1 — 
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dv 2 + 2cdrdv + r 2 dfll 



n—2i 
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where c = ±1 and dCl 2 l _ 2 stands for the metric of the 
unity (n — 2)-dimcnsional sphere, assumed here to be 
spanned by the angular coordinates (81,62, ... , n -2), 
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For the case of an ingoing radial flow, c = 1 and m(v) 
is a monotone increasing mass function in the advanced 
time v, while c = — 1 corresponds to an outgoing radial 
flow, with m(v) being in this case a monotone decreasing 
mass function in the retarded time v. The 4-dimensional 
Vaidya metric with a cosmological constant in the radi- 
ation coordinates has been considered previously in [2l| . 

It is well known that the radiation coordinates are de- 
fective at the horizon [22|, implying that the Vaidya met- 
ric 0]) i s not geodesically complete in any dimension. 
(See [23[ for a discussion about possible analytical ex- 
tensions). Besides, the cross term drdv can introduce 
unnecessary oddities in the hyperbolic equations govern- 
ing the evolution of physical fields on spacetimes with 
the metric ([TJ). Typically, the double- null coordinates 
are far more convenient. This was the main motiva- 
tion of Waugh and Lake's work[24|], where the problem 
of casting the 4-dimensional Vaidya metric with A = in 
double-null coordinates is considered. As all previous at- 
tempts to construct a general transformation from radia- 
tion to double-null coordinates have failed, they followed 
Synge(25| and considered Einstein's equation with spher- 
ical symmetry in double-null coordinates ab initio. The 
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resulting equations, however, have revealed to be not an- 
alytical solvable for generic mass functions. Waugh and 
Lake's work was recently revisited in 26], where a semi- 
analytical approach allowing for generic mass functions 
was proposed. Such approach consists in a qualitative 
study of the null-geodesics, allowing the description of 
light-cones and revealing many features of the underly- 
ing causal structure. It can be used also for more quan- 
titative analyses, indeed, it has already enhanced con- 
siderably the accuracy of the quasinormal modes anal- 
ysis of 4-dimensional varying mass black holes [3, EH] j 
and it can be also applied to the study of gravitational 
collapse [26[. We notice that another method to construct 
conformal diagrams based on a systematic study of the 
null-geodesics was also recently proposed in [27j |. 

Here, we extend the approach proposed in [26| to the 
n-dimensional (n > 2) case and to allow the inclusion 
of a cosmological constant A. Some new exact solutions 
are also presented. For n — 3, notably, a crucial non- 
linearity can be circumvented and the problem can be 
reduced to the solution of a second order linear ordinary 
differential equation, opening the possibility for analyti- 
cal study of large classes of mass functions. We generalize 
one of Waugh and Lake's exact solutions by introducing 
the case corresponding to m(v) oc u"~ 3 , n > 3, describing 
a naked or shell focusing singularity in a n-dimensional 
spacetime. We present also some explicit examples of the 
use of the semi-analytical approach to the study of the 
causal structure of some particular solutions. 

In the next Section, the main equations are derived 
and the exact solutions are presented. Section IIIII is de- 
voted to the introduction of the semi-analytical approach. 
Some explicit examples are presented, particularly the 
two apparent horizons case of an evaporating black hole 
in a de Sitter spacetime and the BTZ black hole under- 
going a finite time interval mass accretion process. The 
last Section is left to some concluding remarks. This 
work has two further appendices. The first one presents 
explicitly the n-dimensional geometrical quantities nec- 
essary to reproduce the equations of the Section [Til an d 
the last contains the polynomial manipulations involved 
in the calculations of the apparent horizons of Sections 
[TT]and[TTTJ 



II. THE METRIC 

The n-dimensional spherically symmetric line element 
in double-null coordinates (u, v, 0\, . . . , n -%) is given by 



dial null-fluid in the eikonal approximation is given by 



Tab = —h(u,v)k a k b , 
an 



(4) 



where k a is a radial null vector. We will consider here, 
without loss of generality, the case of a flow along the 
u-direction. The case of simultaneous ingoing and out- 
going flows in 4 dimensions was already been considered 
in [281 ] . Einstein's equations are less constrained in such 
case, allowing the construction of some exact similarity 
solutions, which, incidentally, can also be generalized to 
the n-dimensional case in the light of the present work. 
Einstein's equations with cosmological constant A 



Rab ~ ^9abR = -J^9ab + ^T ab 



(5) 



implies that, for the for the energy- momentum tensor (|4|). 

In 

R = - *-A. (6) 

n — 2 

Using (|A6p and © , Einstein's equations for the metric 
(|3"|) and the energy-momentum tensor ^ read 
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For n^3, differentiating Eq. (fit))) with respect to u and 
then inserting Eq. leads to 



r n ~ 2 r 12 



2A 



/ (n-2)(n-l) 



= -A, 



(11) 



where A(y) is an arbitrary integration function. The n = 
3 case will be considered below. Now, differentiating Eq. 
(flTJ)) with respect to v and using ([8]) and (fTTj) gives 



h 



n-2\ fA 2 



n — 3 J r n 2 r\ 
Eq. |7]) is ready to be integrated 

/ = 2Bn, 



(12) 



(13) 



where B(v) is another arbitrary integration (and nonva- 
nishing) function. From (fT2")) and (p~5|) . one has 



ds 2 = —2f(u, v)du dv + r 2 (u, v)dQ 2 l _ 2 



(3) 



where f(u, v) and r(u, v) are non vanishing smooth func- 
tions. We adopt here the same conventions of [24( and 
[26] ]. In particular, the indices 1 and 2 stand for the dif- 
ferentiation with respect to u and v, respectively. 

The energy-momentum tensor of a unidirectional ra- 



.= -21^^ 



n — 3 



(14) 



Finally, by using (10]) and (l3"l). Eq. flTTJ) can be written 

as 



t 2 = —B I 1 — 



2A 



2A 



(n-3)r»- 3 (n-2)(ra-l) 



(15) 
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Note that © follows from Eqs. |Q]§ and (15]). Einstein's 
equations are, therefore, equivalent to the equations 
(fT4| . and (fl"5|) . generalizing the results of [24[ an d |2£ 

In order to interpret physically the arbitrary integra- 
tion functions A(v) and B(v), let us transform from the 
double-null coordinates back to the radiation coordinates 
by means of the coordinate change (u, v) — > (r(u,v),v). 
Such a transformation casts ([3]) in the form 



ds 2 = ABrydv 2 - ABdrdv 



(16) 



where (|13|) was explicitly used. Comparing (1} and (fl6|) 
and taking into account (|15|) . it is clear that with the 
choice 



B 



1 A 2 

'2|A 2 | 



(17) 



for ^2 7^ 0, the function A(v ) should represent the mass 
of the n-dimensional solution, suggesting the following n- 
dimensional generalization of the mass definition based in 
the angular com pon ents of the Riemann tensor proposed 
for for A = in [H : 



-r n - 3 R a 



(18) 



k > 1, see (|A5|) and (fT5)) . Note that for constant A, the 
choice of the function B(v) is irrelevant since it can be 
absorbed by a redefinition of v. Note also that the weak 
energy condition applied for ^ requires, from (|14p . that 
BAi < 0. Since B must be nonvanishing from (|13p . A(v) 
must be a monotone function. 

The Kretschmann scalar K — R a bcdR abcd could be also 
invoked to interpret the integration function A. Taking 
into account Eqs. (fT3|) and (fl~5|) . we have from (|A8|) 



case of a charged radial null fluid). In j26[, another 4- 
dimensional exact solution corresponding to A{v) = k/v 
and A = was also presented and discussed. These are 
the only varying mass analytical solutions obtained in 
double-null coordinates so far. We notice, however, that 
Kuroda was able to construct a transformation from radi- 
ation to double-null coordinates for some other particular 
mass functions in four dimensions [lOj. 

In the following Section, we will present a semi- 
analytical procedure to attack the problem of solving 
Eqs. (fl"3")) - (fl"5j) for general mass functions obeying the 
weak e nerg y condition, generalizing in this way the re- 
sults of [2g| obtained for n = 4 and A = 0. The approach 
allows us to construct qualitatively conformal diagrams, 
identifying horizons and singularities, and also to eval- 
uate specific geometric quantities. Before, however, we 
notice that the Waugh and Lake's linear 4-dimensional 
solution can be also generalized for n-dimensions. Let us 
consider the mass function 



A(v) = Xv 



71 — 3 



(20) 



V > 0, A > 0. For this mass function and A = 0, with 
the choice (fTTj). Eq. (JTSJ) reads 



1 



A 



r-2 
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V r 



2 n-3 \r 
which can be integrated as 



n— 3 



(21) 
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2 n-3 



-(n-3) 



ds = D(u), 



where D(u) is an arbitrary integration function, imply- 
ing, from Eq. (fl"3|) . that 



K 



(n-l){n-2) 2 A 2 



(n-3) r z(n~i) ( n _l)( ?1 -2) : 



:A 2 . 



(19) 

It is clear from (fT^|) that r = is a singularity for A ^ 
and that A acts as a gravitational source placed at r = 0. 

The problem now may be stated in the same way of 
the 4-dimensional A = case[26[: given the mass func- 
tion A(v) and the constant B, one needs to solve Eq. 
(fT5|) . giving rise to the function r(u, v). Then, f(u, v) and 
h(u,v) are calculated from p3[) and (|14[) . The arbitrary 
function of u appearing in the integration of (|15p must 
be chosen properly [24| in order to have a non-vanishing 
f(u,v) function from (fl"3"|) . Unfortunately, as stressed 
early by Waugh and Lake[24j. such a procedure is not 
analytically solvable in general. They, nevertheless, were 
able to find some regular solutions for Eqs. (|13p - (fT5|) for 
n = 4 and A = 0, namely the linear (A(v) — Xcv) and a 
certain exponential (A(v) — 4 (aexp(/3cw/2) + 1)) mass 
functions (A, a, and (3 are positive constants, c = ±1, cor- 
responding to ingoing/outgoing flow, respectively). The 
4-dimensional linear mass case was also considered in 
Q in great detail and in a more general situation (the 



/ = 



v n - 2 D 1 



y.71 — 3 



r\«-2 1 /r\»-3 A 
v J 2 \vJ n — 3 



(23) 



The integration function D must be chosen in order to 
have a regular /. As in the 4-dimensional case originally 
considered by Waugh and Lake, we have three qualita- 
tively distinct cases: < A < A£, A = A£, and A > A°, 
where 



AS — 



2(n- 2) 



(24) 



For A > A£, for instance, the quantity inside the square 
brackets in (f2"3"| does not vanish, see Appendix [B] The 
integral in ((22)) can be (numerically) evaluated and, for 
instance with the Waugh and Lake's choice D(u) = —u, 
r(u, v) can be also (numerically) determined. The result- 
ing causal structures are the same ones as the Waugh 
and Lake's 4-dimensional cases. We will construct the 
relevant conformal diagrams for the three qualitatively 
distinct cases as an application of the semi-analytical ap- 
proach. 
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A. The three dimensional case 

In lower dimensional (n < 4) spacetimes, the Riemann 
tensor is completely determined by its traces, namely 
the Ricci tensor and the scalar curvature, implying some 
qualitative distinct behavior for the solutions of Ein- 
stein's equations in these situations [3(J. For n = 3, Eq. 
(fTUll reads 



rr 12 
f 



+ Ar 2 = 0. 



(25) 



It is interesting to compare this last equation with (jTTJ) . 
Eq. (fl"3"|) is still valid and by using it, Eq. (|23|) can be 

easily integrated, 



r 2 + ABr = C, 



(26) 



where C is an arbitrary integration function which we 
call, by convenience, C(v) — —B(v)A(v), leading to 



r 2 = —B (-A - Ar 2 ) , 



(27) 



which corresponds to the 3-dimensional counterpart of 
Eq. ([T5|). Eq. © for n = 3 also follows from ||T3J) and 
(|2T|) . As for the n-dimensional case, by using (fHij) and 
(|27|) . one gets from (J8]) the last equation 



h = -B^. 

r 



(28) 



One can show that the integration functions A and B 
have the same physical interpretation as the higher di- 
mensional cases by introducing the radiation coordinates. 
We have for the 3-dimensional case 



ds 2 



-(2B) 2 (-A-Ar 



dv 1 



4Bdrdv + r 2 d6 2 . (29) 



It is clear that with the choice (fTT|) for B, the function 
A plays the role of the BTZ black- hole mass[3l| (A < 
in this case). Note, however, that for n = 3 there is no 
purely angular components of the Riemann tensor and, 
therefore, there is no equivalent of the mass definition 
(I18p . The Kretschmann scalar in this case reads simply 



K = 12A 2 . 



(30) 



In fact, if terms involving distributions are take into ac- 
count properly, curvature invariants as K reveal that the 
origin for the BTZ black hole is a conical singularity (See, 
for a recent rigorous analysis, [32|). 

In contrast with the n-dimensional case, the non- 
linearity present in (|27p can be easily circumvented. Let 
us consider A 2 > 0, B = -1/2, and A = -I /I 2 (Other 
cases follow straightforwardly) . By introducing the linear 
second order ordinary differential equation 



it is easy to show that 
r(u,v) = -2£ : 



2 P(U)<(V) 

P(u)w a (v) + w b (v) ' 



(31) 



(32) 



is solution of (|27|) . where w a (v) and Wb(v) are the two 
linearly independent solutions of (|3T|) . and P(u) is an 
arbitrary function. The second order linear equation (|3Tj) 
has analytical solution in closed form for many functions 
A(v). For instance, for A(v) oc v a (or A(v) oc exjpav), 
a 6 E, Eq. (|31|) is equivalent to the Bessel equation, 
after an appropriate redefinition of v (v — > ti( Q / 2+1 ) or 
v — > e au / 2 ). Even for the case where no solution in closed 
form can be obtained, the main analytical properties of 
w(v) can be inferred easily from pip due to its linearity. 



III. SEMI-ANALYTICAL APPROACH 

We review here the semi-analytical approach proposed 
in [2(| for the solution of the equation (fT5|) (or, anal- 
ogously, l|27p for the 3-dimensional case). Firstly, no- 
tice that in double-null coordinates the light-cones corre- 
spond to the hypersurfaces with constant u or constant 
v. One can, in this case, deduce the causal structure of a 
given spacetime by considering the set of null geodesies 
corresponding to u = constant or v — constant. The 
function r(u, v), obtained as the solution of (|15p with a 
careful choice of initial conditions [241. [26j . has also a clear 
geometrical interpretation: it is the radius of the (n — 2) 
dimensional sphere defined by the intersection between 
the hypersurfaces u — constant and v — constant. Once 
r(u,v) is obtained from (|15p . f(u,v) and h(u,v) can be 
calculated directly from Eqs. (fTB")) and (fl"4")) . giving all 
the information about the spacetime in question. Inci- 
dentally, Eqs. (fT5")) along constant u, and, consequently, 
along a portion of the light-cone, is a first order ordi- 
nary differential equation in v. One can evaluate the 
function r[u, v) in any spacetime point by solving the v- 
initial value problem knowing r[u, 0). For instance, the 
usual constant curvature empty spacetimes (^4 = 0) can 
be obtained by choosing r(u, 0) = u/2, leading to 



2A 



(n - 2)(n - 1) 



ds = -(v + u). (33) 



As it is expected, for the n-dimensional de Sitter case 
(A > 0), r = y/(n — 2)(n — 1)/2A corresponds to a cos- 
mological horizon. The condition ri (u. 0) ^ is sufficient 
to assure that f(u,v) ^ everywhere [26] . 

Let us focus now on the solutions r(v) — r(u, v) of (|15p . 
considered ELS Si Cl first order ordinary differential equation 
in v for constant u — u, keeping in mind that they indeed 
describe how r(u, v) varies along the u — u portion of 
the light-cone and, thus, that they are closely related to 
the causal structure of the underlying spacetime. This 
situation is shown schematically in Fig. [T] 

For n > 3, the curves defined by 



2A 



n— 3 



2A{v) 



(n-2)(n- 1) 



(34) 



correspond (if they indeed exist) to the frontiers of cer- 
tain regions of the (v, r) plane where the solutions r{v) 



FIG. 1: The solutions r(u,v) of (|15[) for the Schwarzschild 
case (A(v) = m, see H3j). The dotted lines correspond to 
the lines of constant r. In the region I (the interior region), 
the values of r decreases monotonically upward, varying from 
r — 2m, the horizon corresponding to the degenerated hy- 
perbole u = Uo and v = vo, to the singularity r = 0. In 
the exterior region (II), on the other hand, the values of r 
increases monotonically leftward, varying from r = 2m to the 
spatial infinity r — ex. Some slices of constant u = u are 
depicted by the lines with arrows. The values of r(u, v) read 
along such lines correspond to the solutions of (|15[) with a 
certain initial condition r(0) = r«. Since the lines of constant 
u are portions of the light cone, the function r(u, v) describe 
how the light cones cross the (n — 2)-dimensional spheres of 
radius r. In the interior region, for instance, all null geodesies 
eventually reach the singularity, whereas in the exterior, the 
null geodesies with constant u escape to the null infinity. Be- 
tween such regions, there is an event horizon at u = uq. We 
notice that all this information is coded in the qualitative 
behavior of the solutions r(u,v) for the first order ordinary 
differential equation corresponding to (|15[) considered along 
constant u — u. This is the basic idea of our semi-analytical 
approach. 



of (I15p have qualitative distinct behaviors. For the de 
Sitter case, can define two non intersecting curves, 
whereas for the anti de Sitter case (A < 0) there is only 
one curve. (See Appendix |B|) . Let us suppose, for in- 
stance, a de Sitter case with A 2 (v) > (and B = —1/2. 
Outgoing radiation flows, 3-dimensional spacetimes, or 
the anti de Sitter case follow in a straightforward man- 
ner). Let us call rEn(v) and rcn(v) the curves defined 
by ([3i)l. see Fig. [2] Due to ([it]), for all points of the 
plane (v,r) below rEH(v), r 2 < 0. Hence, any solution 
r(v) entering in this region will, unavoidably, reach the 
singularity at r = 0, with finite v. Suppose a given so- 
lution rj(u) with initial condition r;(0) = n enters into 
the region below teh{v)- As, for smooth A, the unique- 
ness of solutions for (|15[) is guaranteed for any point with 
r > 0, any solution starting at r(0) < r; is confined the 



FIG. 2: In the region below the curve rEH(v) (the apparent 
event horizon), T2 < 0. Any solution r(v) of (|15[) that enters 
into this region will reach the singularity at r = with fi- 
nite v. Solutions confined to the ri > region always escape 
from the singularity and tend asymptotically to rcniv) (the 
apparent cosmological horizon). The region above rcn(v) 
corresponds to the cosmological exterior, any solution there 
also tends asymptotically to rcn(v). The dashed line would 
correspond to the apparent event horizon for the case of a 
vanishing A (The unique solution of (|34|l for A = 0). The 
conforms! diagram is inserted. (The case depicted here cor- 
responds to: n = 5, A = 1, and A(v) = (6/5)tanhn, v > 0. 
In this work, we do not consider maximal extensions. For a 
discussion of maximal conformal diagram for de Sitter and 
anti-de Sitter black holes, see [jj^], for instance. ) 

the region below ri(v) and will also reach the singularity 
at r = 0, with finite v. On the other hand, suppose that 
a given solution r c (v) with initial condition r c (0) = r c 
never enters into the region bellow the curve teh(v). So- 
lution starting at r(0) > r e are, therefore, confined to the 
region above r c (v) and will escape from the singularity 
at r = 0. Hence, one has two qualitatively distinct be- 
havior for the light-cones, see Fig. [5J For r(0) < r- u the 
future direction points toward the singularity and all null 
geodesies eventually reach r = 0. This is the typical situ- 
ation in the interior region of a black-hole. For r(0) > r c , 
the u = constant portion of the light cones escape from 
the singularity. The curve teh{v) plays the role of an 
apparent event horizon. For a given v, all solutions r(v) 
in Fig. [2] such that r(v) < teh{v) will be captured by 
the singularity. Solutions for which r(v) > rEii(v) are 
temporally free, but they may find themselves trapped 
later if Teh{v) increases. The event horizon corresponds 
to the last of these solutions trapped by the singularity 
and is located somewhere between Ti{v) and r e (v). 

In the region between rEti{v) and tch{v), the solu- 
tions r(v) have r 2 > 0. They escape from teh{v) and 
tend asymptotically to rcn(v), which plays the role of an 
apparent cosmological horizon. The region above rcH{v) 
corresponds to the cosmological exterior (See, for in- 
stance, [34|. For ^4 = 0, the cosmological exterior region 
corresponds to the choice ro = oo in (|55|) ). There, r 2 < 
and the solutions also tend asymptotically to tch{v)- 
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The case where teh{v) and rcn(v) converge to the same 
function gives rise to the interesting n-dimensional Nariai 
solution [35J. 

Summarizing, since the partial differential equation 
(|15p for r(u, v) does not involve u explicitly, it is always 
possible to discover how the light-cones cross the (n— 2)- 
dimensional spheres of radius r by exploring, in the plane 
(v,r), the solutions r{u,v) of the first order differential 
equation corresponding to (|15[) considered on the con- 
stant u null-geodesics. From these informations, one can 
reconstruct the causal structure of the underlying space- 
time. For typical cases, we will have some spacetime 
regions where the light cones necessarily end in the sin- 
gularity at r = 0, whereas for some other regions it is 
always possible to escape to the null infinity. The fron- 
tier of these regions must correspond to event horizons. 
Let us illustrate the approach with some explicit exam- 
ples. 

A. Waugh and Lake generalized solutions 

As the first application of the semi-analytical ap- 
proach, let us consider the generalized Waugh and Lake 
solutions corresponding to the choice ([20]) for the mass 
function, n > 3 , and A = 0. From the discussion of 
Section |TT1 we expect three qualitatively distinct case ac- 
cording to the value of A. In fact, these solutions have 
the same causal structure for any n > 3. The frontier 
of the region in the (v, r) plane where all the solutions 
of (pnj) reach the singularity at r = corresponds to the 
straight line 

ro(v)= f-^y~'v. (35) 
Taking the ^-derivative of (|2ip. one gets 

r22 = A ^(- 1+ 2r^3(r) J' (36) 

The regions in the plane (v, r) where the solutions obeys 
r22 = are the straight lines defined by 

(-) -i(-) +-^ = 0- (37) 
\v/ 2 \vJ n — 3 

One can now repeat the analysis done for the n = 4 
case by considering the three qualitative different cases 
according to the value of A and the possible solutions of 
(|3"T|) . For this purpose, we will consider the solutions of 
(121 p with the initial condition r(u, 0) oc u. 

For A > A£, with \ c n given by Eq. J37]) has no 

solution and T2i < for all points with v > 0. The only 
relevant frontier in the plane (v,r) is the ro(u) straight 
line. All solutions of (f2Tj) are concave functions and cross 
the r-i — line, reaching the singularity with a finite v 
which increases monotonically with u. The causal struc- 
ture of the corresponding spacetime is very simple. There 



is no horizon, and all future cones end in the singularity 
at r = 0. The associated conformal diagram is depicted 
in Fig. [3] 




FIG. 3: Conformal diagrams corresponding to the mass func- 
tion (|20p and A = 0. Cases (b) and (c) exhibit naked and shell 
focusing singularities. For all cases, the curves in the (v,r) 
are similar to those ones presented in [2(|, corresponding to 
the n — 4 case. 

For A = A£, r-22 = only on the line 

Mv) - £^v. (38) 

This straight line itself is a solution of pTjl . All other 
solutions are concave functions. Note that, in this case, 
Th{v) > i"o(v) for v > 0. We have two distinct qualitative 
behavior for the null trajectories along u constant. All 
solutions starting at r(0) > arc confined to the region 
above rh(v). They never reach the singularity, all tra- 
jectories reach X + . However, in the region below r^(u), 
we have infinitely many concave trajectories starting and 
ending in the (shell-focusing) singularity. They start at 
r(0) = 0, increase in the region between r^(u) and r (v), 
cross the last line and reach unavoidably r — again, 
with finite v. The trajectory rh(v) plays the role of an 
event horizon, separating two regions with distinct qual- 
itative behavior: one where constant u null trajectories 
reach I + and another where they start and end in the 
singularity. This behavior is only possible, of course, be- 
cause the solutions of (f!?T]) fail to be unique at r = 0. The 
relevant conformal diagram is also shown in Fig. [3] 

For < A < A° , we have three distinct regions accord- 
ing to the concavity of the solutions. They arc limited by 
the two straight lines r_(u) and r+(v) defined by (|37[) . 
We have ro(v) < r_(«) < r + (v) for v > 0. Both straight 
lines r_(u) and r + (v) are also solutions of (JUJ). Between 
them, solutions are convex. Above r + and below r_ , so- 
lutions are concave. The last line is the inner horizon. 
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Inside, the null trajectories with constant u start and 
end in the singularity r — 0. The line r + is an outer 
horizon. Between them, the constant u null trajectories 
starts in the naked singularity and reach I + . Beyond the 
outer horizon, all constant u null trajectories escape the 
singularity, see Fig. (J3j> . 



< u < 0. For u > uq, we have a empty n-dimensional 
de Sitter spacetime. We can now apply the same proce- 
dure we have used for the previous cases. We call teh{u) 
and rcniu) the curves on the (u,r) plane delimiting the 
regions where the solutions of ([42]) have qualitative dis- 
tinct behavior. (See Fig. 2]). In this case, since A% < 0, 



B. Black- hole evaporation in a de Sitter spacetime 

As an example of using the semi-analytical approach 
for A ^ 0, let us consider the case of a evaporating 
black hole due to Hawking radiation in a de Sitter n- 
dimensional (n > 3) spacetime. The mass decreasing 
rate for these n-dimensional black holes can be obtained 
from the n-dimensional Stefan-Boltzmann law (see, for 
instance, [11]), leading to 



dm 
du 



-am 



-2/(n-3) 



(39) 



where a > is the effective n-dimensional Stefan- 
Boltzmann constant. Eq. (|39|) can be immediately in- 
tegrated, leading to 



m(u) = mo 1 



Uq 



(40) 



where uq in the lifetime of a black hole with initial mass 
m , 



auo 



n - 3 (n-l)/(n-3) 



n - 1 



(41) 



We recall that (f39|) is not expected to be valid in the 
very final stages of the black hole evaporation, where 
the appearance of new emission channels for Hawking 
radiation (37l| can induce changes in the value of the con- 
stant a, and maybe even the usual adiabatic derivation 
of Hawking radiations is not valid any longer. We do not 
address these points here. We assume that the black hole 
evaporates completely (|40[) for all < u < uo and that 
m = for u > uq. Notice that, since m' — > 00 for u — ► uq, 
a naked singularity will be formed after the vanishing of 
the black hole [3 . 

We can consider now a n-dimensional Vaidya met- 
ric with a mass term of the form given by (|40jl . This 
is the first case of an outgoing radial flux we consider, 
and according to (|T7|) . B = 1/2. However, in this case, 
f(u, v) has different sign and the temporal directions for 
the B = —1/2 case are now spacelike. The transforma- 
tion (u, v) — > (y, —u) restores the temporal and spatial 
directions [11| . One needs to keep in mind that, for B < 
(ingoing radiation fields), v corresponds to the advanced 
time, whereas for B > (outgoing radiation fields), u 
should correspond to the retarded time. In this case, the 
equivalent of Eq. ([To]) for the mass function (j40|) reads 



3 1 



2mo (1 — u/uqY 



2A 



(n- l)(n-2) 




FIG. 4: The (u, r) plane for an evaporating black hole in a 
n-dimensional de Sitter spacetime. The apparent horizons 
teh(u) and rcniu) delimit the regions where the solutions 
of (|42p have qualitatively distinct behavior. Since the black 
hole vanishes for u = uq, for u > uo the only remaining 
horizon corresponds to tch(u). The conformal diagram is 
inserted. For u < uo, it represents an usual black hole. For 
u > uo, we have the empty space left after the vanishing of 
the black hole. The dashed line HC (extension of the black 
hole horizon EH) is a Cauchy horizon. Beyond it, we have a 
breakdown of predictability due to the naked singularity left 
by the black hole. Part of the cosmological exterior region (the 
region above CH and below HC) is completely insensitive 
to the singularity. (The case depicted here corresponds to: 
n — 6, A = 5/2, and mo = uo = 1, u > 0.) 



the curve Teh (v) decreases and reach r = for u = uq, 
corresponding to the total vanishing of the black hole. 
After this point, tch{v) is the only horizon (a cosmo- 
logical one) of this spacetime. Before u = uq, however, 
some solutions of (l4"2"j) that crossed teh{v) could reach 
the singularity at r = 0. The exterior cosmological region 
(region of the (it, r) plane above tch(v)) is insensitive to 
the singularity. The resulting conformal diagram shown 
in Fig. 2] It represents a spacetime with a cosmological 
horizon and, in its interior, a black hole that evaporates 
completely leaving behind a naked singularity. 



(42) 



C. Mass accretion by BTZ black-holes 

As an explicit example of the n = 3 case, let us consider 
a BTZ black hole[3l| undergoing a mass accretion process 
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during a finite time interval. The C -class mass function 



mi, V < —Vq, 

A(v) — ^ av (3«q — v 2 ) + b, —v < v < v , 

Hlf, V > Vq. 



(43) 



does correspond to such a situation. Here, m\ and TOf 
stand for, respectively, the initial and final mass, and the 
constants a — (mt — mj)/4wg and b = (n%f + m{)/2 are 
chosen in order to have smooth matches at v = ±«o. 
For the BTZ case, A = — 1/£ 2 and the apparent horizon 
corresponding to the region of the (r, v) plane where r<i = 
is given simply by t"eh(^) = £\/A(v), see (|2"T)) . Since 
A 2 > 0, we adopt B = -1/2. The solutions of J27|) below 
the curve teh(^) have r 2 < 0, whereas those ones above 
have T2 > 0, see Fig. [5] For u < — v , and analogously 




FIG. 5: The (v,r) plane for a BTZ black hole accreting mass 
according to ()43|) . In the region below the curve rEH(v) (the 
apparent event horizon), T2 < 0. The constant u = u null 
geodesies associated to the solutions r(u, v) of (|27|) entering 
into this region reach the conical singularity at r — with 
finite v. The null geodesies confined to the r2 > region can 
escape from the singularity, reaching the null and spatial in- 
finity. The conformal diagram is inserted. (The case depicted 
here corresponds to: A = — 1, B = —1/2, and mi = mf/9 = 1. 

) 



for v > vq, Eq. (|2"T| can be easily integrated 



r{u, v) — lyfmi 



r(u, v) + l^/mi 



Di(u)e 



/m,iV I i 



(44) 



where Di(u) is a monotone arbitrary integration function. 
The region r > t^/mi, is obtained by choosing Di{u) > 0, 
whereas r < l^/m{ corresponds to D\(u) < 0. In the inte- 
rior region of a BTZ black hole, we see from (14"4"1) that the 
constant u null-geodesics reach the conical singularity at 
r = in a finite time given by (l/£) ln(— 1/Di(u)). In con- 
trast with the n > 3 cases, the solutions r(v) here do not 
reach r — perpendicularly (See Fig. [5]). For the exterior 
region, the null-geodesics can reach the spatial infinity 
r = oo also in a finite time given by (l/£) ln(l/Dj(u)). 

The solutions r(v) for the mass function (|4"3"|) are ob- 
tained by matching, at v = ±vq, the BTZ black hole (14~4"1) 



with the accretion process corresponding to — vq < v < 
vq. This situation is depicted in the Fig. [5] 



IV. CONCLUSION 

We considered here the problem of constructing the Tri- 
dimensional Vaidya metric with a cosmological constant 
A in double- null coordinates, generalizing the work [26| . 
where a semi-analytical approach is proposed to attack 
the equations derived previously by Waugh and Lake 24 
for n = 4 and A = 0. Some new exact solutions are 
also presented. For n — 3, in particular, the problem 
reduces to the solution of the second order linear ordinary 
differential equation ([3T|) . allowing the analytical study 
of large classes of mass functions A(v). As an example, 
for A < and A(v) — mo + mitanhu, too > mi, Eq. 
(|3ip is equivalent to the hypergeometric equation, after 
a proper redefinition of v. This situation corresponds to 
a BTZ black hole with mass too — mi that accretes some 
mass smoothly and ends with mass too + mi. This kind 
of smooth exact solution, with two "asymptotic BTZ" 
regions (v — * ±oo) is very interesting to the study of 
creation of particles in non stationary spacetimes [38j . 

The generalized semi-analytical approach is also use- 
ful to the study of quasinormal modes of varying mass 
?i-dimensional black holes. In particular, an interesting 
point is the study of the stationary regime for quasi- 
normal modes, described in [14], in higher dimensional 
space and in the presence of a cosmological constant. The 
semi-analytical approach can help the understanding of 
the relation anticipated by Konoplya 39] for the quasi- 
normal modes of n-dimensional black holes with A = 0: 
lur oc Tir^ 1 , where ton and ro stand for the oscillation fre- 
quency and the black hole radius, respectively. Some pre- 
liminary results [lH suggests the existence of a stationary 
regime such that, for slowly varying mass n-dimensional 
black holes, w_r(w) oc tiTq (i>). As another application 
for A ^ 0, we have the dynamical formation of a Nariai 
solution, which could be accomplished by choosing A(v) 
such that teh{v) — > tch{v) for larges v. These points 
are still under investigation. 
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APPENDIX A: SOME GEOMETRICAL 
QUANTITIES 

We present here some geometrical quantities necessary 
to reproduce the equations of Section II. Such quantities 
are the n-dimensional (n > 2) generalization of those 
ones calculated previously by Synge[25[ and Waugh and 
Lake[24] for n = 4. Only nonvanishing terms are given. 



9 



Recall that our n-dimensional line element is given by 



ds 2 = —2f(u, v)du dv + r 2 (u, v)dfl 2 l _ 2 



(Al) 



where dfl 2 l _ 2 is the unity (n — 2) sphere, assumed here to 
be spanned by the angular coordinates 6i, i = 1 . . . (n — 
2), 



Re* 



IL 

\l=l 



fc-1 

n s[n20t 



1 + 2 



Using (|13j) . we get from the purely angular components 
of the Riemann tensor: 



n-2 I i-l 



i=i \j=i 

The Christoffel symbols are 



f 

e k _ n 



1iU 



vv ~ 


/a 
" /' 


P#fc 


_ ^2 

r 


fc-1 




»=i 




fc-i 




J] sin 2 

f=l 




/ fc-i 






sin 2 6>i 







(A2) 



(A3) 



The totally covariant Riemann tensor R a bcd is given by 
R _ /1/2 _ , 



'fc-1 



Ru9 k u9 k 
Rv9 k v9 k 



T \ 1 1 Sm 

\t=l 

n >>n ' 

i=i 

/fc-i 



/ 

/ 



?"22 



Ru9 k e k v = r II sin2 ^ r 12 



''2 



R 9 1 9 k 9 1 ~ 1+ B> 



(A5) 



for any k > 1. As for the n = 4 case[29[, such components 
are invariant under transformations involving solely the 
double-null coordinates (u,v) — > (f/(u,u), V(tt,u)) and 
are candidates for defining the mass of the solution for 
A = 0. 

The Ricci tensor, obtained by the contraction R a b = 

R C ab is s iven b y 



Ru 
R u 
R,, 



(n-2) 



fin _ ni 

fr r 



(A6) 



/1/2 /12 , <^ r 12 
-(n-2) — 



/ 2 / 
(n-2) 
/fc-i 



■Re 



n 

vi=l 



"2 



/27"2 r22 

r 

(rri2 + (n — 3)rira) + (n — 3) 



/ 



(A4) 



Two scalar quantities are relevant for our purposes: 
the scalar curvature R = g ab R a b and the Kretschmann 
scalar K — R a bcdR ahcd ■ They are given by 



R 



and 



/1/2 



2/ -i^ + (n-3)(l + ^ 2 
(A7) 



K 



(n-l)(n-3) 



4 
P 



72 r 12 + r ll r 22 + 



nr-if\h r-nnf\ ni^/a 



1'2 



^1^2 



4 

7* 



/ 2 



/ 2 

2 , ,2 
' ./l2 



- 2 



/ 

/1 2/1/2 
/ 



/ 



(A8) 



APPENDIX B: SOME POLYNOMIAL RESULTS 

The evaluation of the zeros of / given by (|2"3"|) involves 
finding the positive roots of the polynomial 

p n {x)=x n - 3 [x-^j =-^3, (Bl) 



for A > and n > 3. The only roots of p n (x) are x = 1/2 
and ie = 0, the latter with multiplicity n — 3. For posi- 
tive IE, p n (x) < only in the interval x E (0, 1/2). The 
minimal value of p n {x) in such an interval corresponds 
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to 



-A' 



1 / n-3 



n-3 n n-3 \2(n - 2) 



where 



c p = 

max 



n — 3 
2(n-2) : 



(B2) 
(B3) 



the unique root of p' n (x) in the interval. Thus, for A > A^, 
the polynomial (|Bip has no roots. For A — A^ there is 



only one root (x 



c ), while for < A < A^j, there are 



two roots xi and X2, with < Xi < x^ ax < X2 < 1/2. 

The determination of the curves ((34)) involves the eval- 
uation of the positive roots of the polynomial 



Qn(x) 



n — 3 



2A 



(n- l)(n-2) 



2.4 



(B4) 



with A > and n > 3. For 
the only root of q n {x) is x 
positive for all x > 0, and 
positive root. Equation 
one curve. 

On the other hand, for 
q n (x) has always 3 roots: x 
x = 0, the last with multipl 
itive x, q n (x) > only for 



the anti de Sitter case, A < 0, 
= 0, q n (x) and q' n {x) are both 
therefore, (|B4[) has only one 
in this case, defines only 

the de Sitter case, A > 0, 
= ±y/(n- l)(n- 2)/2A and 
icityn — 3. Moreover, for pos- 
x e [0,V(n-l)(n -2)/2A], 



and the maximum value of q n (x) in such an interval is 

fc(O) where 



(n-2)(n-3) 



2A 



(B5) 



the only point of the interval where q' n (x) = 0. For 
< 2A/(n — 3) < g„(x max ), (|B4|) has always two roots 
xi and x 2 , X\ < x max < x 2 . For 2A/(n - 3) = <7«(x max ), 
there is only one root x = x max , and for 2A/(n — 3) > 
p„(x^ ax ) there is no root at all. In this case, provided 
that < 2A/(n - 3) < g„(x max ), Eq. ([34]) defines 
two curves r EH (v) and r C H( v ) (See Fig. [5]), such that 
teh(v) < x« lax < r C H{v). 

If r(i>) is a curve defined by (|34[) . its derivative is given 
by 



n-3 r 



4Si <»-»>- 



2A 

n- 2 ? 



(B6) 



If A < 0, it is clear that r'(«) has the same sign of A'{v). 
For A > 0, the curves rEH{v) and tch{v) have qualita- 



tively distinct behavior. Since teh{v) < x^ 



EH 



has the same sign of The curve rcii(v) has the 

converse behavior, since tch(v) > x max . 
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